LARGE N c PHYSICS FROM THE LATTICE 



M. TEPER 

Theoretical Physics, University of Oxford, 1 Keble Road, Oxford 0X1 3NP, UK 
E-mail: teper&thphys. ox. ac. uk 

I summarise what lattice methods can contribute to our understanding of the 
phenomenology of QCD at large N c and describe some recent work on the physics 
of SV(Nc) gauge theories. These non-perturbativc calculations show that there 
is indeed a smooth N c — + oo limit and that it is achieved by keeping g 2 N c fixed, 
confirming the usual diagrammatic analysis. The lattice calculations support the 
crucial assumption that the theory remains linearly confining at large N c . Moreover 
we see explicitly that N c = 3 is 'close to' N c = oo for many physical quantities. We 
comment on the fate of topology and the deconfining transition at large N c . We find 
that multiple confining strings are strongly bound. The string tensions of these k- 
strings are close to the M(-theory)QCD-inspired conjecture that cj. oc sin(-7rfc/A r c ) 
as well as to 'Casimir scaling', cr^ oc k(N c — k), with the most accurate recent 
calculations favouring the former. We point out that closed fc-strings provide a 
natural way for non-perturbative effects to introduce 0(1/N C ) corrections into the 
pure gauge theory, in contradiction to the conventional diagrammatic expectation. 

1 Introduction 

As we have seen at this meeting, large-A^ arguments are very useful in 
illuminating many aspects of QCD. Although we may not know very much 
about the detailed physics at N c = oo, we can get a long way by assuming that 
there is a smooth large N c limit which is confining and that N c — 3 is 'close 
to' that limit. Analysis of the colour flow in Feynman diagrams tells us that 
this limit is achieved by keeping constant the 't Hooft coupling, A = g 2 N c , 
and that the leading corrections are 0(1/N C ) in QCD and 0(1/N%) in the 
gauge theory. 

The question I want to address in this talk is: are the above assumptions 
correct and does the colour flow counting survive if we go beyond diagrams 
to a fully non-perturbative calculation? The technique I shall use is to dis- 
cretise the theory onto a space-time lattice, calculate various mass ratios via 
computer simulation, do this for a large enough range of lattice spacings that 
one can confidently extrapolate to the continuum limit; and finally, repeat the 
exercise for a large enough range of N c that one can confidently control the 
approach to N c = oo. 

I will do this for SU (N c ) gauge theories with no quarks. These are cal- 
culations that can be - and have been - done on workstations. The states 
of this theory are purely gluonic so one may call them 'glueballs'. If we take 
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the N c — > oo limit of this glueball spectrum, what we obtain is the glueball 
spectrum of QCD^ c=00 , since we expect no mixing between glueballs and 
quarkonia at leading order in 1/N C (at least for m q ^ 0). The next step 
would be to do QCDn c in the 'quenched' approximation, where all quark 
vacuum bubbles arc neglected. This can be regarded as a relativistic valence 
quark approximation to the theory. At any fixed m q ^ it has QCD^ a=00 
as its N c — > oo limit, since quark vacuum bubbles do not appear at leading 
order in 1/N C , and so it can be used to determine the quarkonium physics of 
that theory. Such a calculation would be interesting and should be possible 
using a Teraflop computer of the kind that is becoming available to a number 
of lattice groups. If however one wants to look at the chiral limit at large N c , 
then one needs to include quark loops and these are calculations that are not 
for the near future. 

In the next section I briefly remind you how one calculates masses us- 
ing lattice simulation and I give you an explicit example to demonstrate that 
such calculations are indeed possible. I then proceed to describe the results 
of a calculation □ of the lightest few glueballs which shows that the approach 
to N c = oo is remarkably rapid: even SU(2) is 'close to' SU(oo) for many 
quantities. These calculations provide some explicit evidence that linear con- 
finement survives at large N c and that the limit is indeed achieved by,keeping 
g 2 N c fixed. A much larger calculation of this kind is now in progress 13. I then 
summarise the lattice results for the tensions of fc-strings and what they 
imply for various model/theoretic expectations. These string tension calcu- 
lations are also interesting because they provide an explicit example of how 
non-perturbative physics may violate the usual large N c diagrammatic colour- 
counting results. Finally I briefly summarise what we are learning about the 
deconfining transition and about topology. 

2 Calculating Masses 

To calculate a mass we construct some operator cf>(t) with the quantum num- 
bers of the state and then use the standard decomposition of the Euclidean 
correlator in terms of energy eigenstates 

c(t) = (^(t)^(o)) = j2\(^\m\ 2 cM-E n t} (i) 

71 

where \n) are the energy eigenstates, with E n the corresponding energies, and 
|f2) is the vacuum state. To be able to evaluate the corresponding Euclidean 
Feynman Path Integrals we discretise continuous space-time to a lattice and 
truncate the infinite volume to a finite hypertorus. We now have a finite 
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number of degrees of freedom and can calculate Feynman Path Integrals using 
standard Monte Carlo techniques. 

The lattice degrees of freedom are SU(A C ) matrices that reside on the 
links of the lattice. In our partition function the fields are weighted with 
exp{S*} where S is the standard plaquette action 



s = -eY,( 1 -w ReTrU p)> ( 2 ) 



N, 

and Up is the ordered product of the matrices on the boundary of the plaquette 
p. For smooth fields this action reduces to the usual continuum action with 
j3 = 2N c /g 2 . However the fields that dominate the Feynman Path Integral 
are rough, all the way to the scale of the lattice spacing a. For these fields we 
can define a running lattice coupling <7z,(a) which reduces in the continuum 
limit to a coupling g(a) in our favourite scheme: 

P 9l(a) ~~ ' 9 2 {a) W 

So by varying the inverse lattice coupling (3 we vary the lattice spacing a. 

If we use a lattice action with reflection postivity, such as the simple 
plaquette action in eqn(|^), then the decomposition in eqn(Q) remains valid, 
except that now t = ant-, so that we obtain the energies from eqn(Q) as aE n 
i.e. in units of the lattice spacing. 

Having calculated some masses ami at a fixed value of a we can remove 
lattice units by taking ratios: ami/ am j — rrii/mj. This ratio differs from the 
desired continuum value by lattice corrections. For our action the functional 
form of the leading correction is known to be 0{a 2 ). Thus for small enough 
a we can extrapolate our calculated mass values 

7-4 = 7t4 + ca z m%(a) (4) 

mj (a) m^O) 

where c depends on i,j and k and the a-dcpcndcnce of m^(a) will make 
differences at (9(a 4 ). At this point we have obtained the mass ratios of the 
continuum theory which is the ultimate goal of our lattice calculations. 

If we want to calculate the lightest mass using eqn(^) then it is clear that 
we have to go to large enough t that the contribution of the excited states has 
died away and the correlation function has acquired a simple exponential fall- 
off with t. At large t, however, the value of the correlation function becomes 
very small and it is not obvious that a numerical approach, with finite errors, 
will be accurate enough. To demonstrate that it can be, I show in Fig|l] the 
correlation function used to extract the lightest SU(4) J PC — ++ glueball 
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Figure 1. SU(4) correlation function for the ++ glueball; exponential fit shown. 



mass in an ongoing calculation 13. On this plot a simple exponential is a 

straight line and it is clear that the corresponding mass can be determined 

very accurately. It is also clear that the simple exponential decay already 

starts at small t. This means that the operator we are using must be a good 

approximation to the lightest glueball wavefunctional. This is no accident; it 

has been obtained by a variational procedure which is an important ingredient 

in the successful lattice calculation of glueball masses, but one which I have 

no time to describe further here □. 

We have just seen that lattice calculations of masses are indeed possible. 

But are they accurate enough to permit a controlled extrapolation to the 

continuum limit? To demonstrate that the answer to this question is yes, I 

Slot in Fig.^Jthe (preliminary) ++ masses obtained from the same calculation 
at four different values of a. The masses have been expressed in terms of the 
confining string tension, a 2 a, which has also been also calculated, and the ratio 
is plotted against a?a. As we note from eqn(Q) the leading lattice correction 
is 0(a 2 ) which means that the continuum extrapolation at sufficiently small 
a will be a straight line - as shown in the plot. This provides an example of 
a typical continuum extrapolation. 
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Figure 2. The lightest SU(4) scalar glueball mass, rrtg + +, expressed in units of the string 
tension, <r, plotted against the latter in lattice units. The a — + continuum extrapolation, 
using a leading lattice correction, is shown. 



3 SU(N C ) Glueball Masses 

In i we calculated the lightest and hrst excited ++ glueball masses and the 
lightest tensor 2 ++ glueball mass. We took the ratio to the string tension 
and extrapolated to the continuum limit as described in Section |^. We did 
this for SU(2), SU(3), SU(4) and SU(5) gauge theories. In FigJ|l plot these 
continuum mass ratios against 1/-/V,?. We expect the leading correction at 
large N c to be 0(l/iV c 2 ), 



rn ; 



N c 



rrij 
rn-, 



NV 



(5) 



which is a simple straight line on our plot. Remarkably, as we see in Fig.|^, the 
mass ratios for all values of N c can be described by just the leading correction 
and the corresponding coefficients are modest in magnitude. 

This shows us that there is indeed a smooth large iV c limit with a finite 
confining string tension. Moreover for these quantities SU(3) is clearly close 
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Figure 3. The mass of the lightest scalar glueball, m ++ , .... expressed in units of the string 
tension, a, is plotted against The N c — > co extrapolation is shown. 



to SU(oo). Indeed, so is SU(2). That is to say, SU(7V C ) gauge theories are 
close to SU(oo) for all values of N c . 



4 't Hooft Coupling 

We have seen that there is a smooth large- N c limit. Is it achieved by keeping 
constant the 't Hooft coupling, A = g 2 N c , as suggested by the standard anal- 
ysis of diagrams EJ? In D=2+l the coupling g 2 has dimensions of mass and 
the question is simply whether g 2 N c /y/a r goes to a non-zero finite constant as 
N c — ► oo. The answer is found to be yes B. Here in D=i+ 1 the coupling runs 
and is dimensionless. The question therefore becomes BQ: is the smooth large 
N c limit achieved by keeping fixed the running 't Hooft coupling, as defined 
on some scale I that is fixed in units of some quantity that partakes of the 
smooth large- N c limit, such as the string tension? To test this we use eqn(||) 
which tells is that a suitable defintion of a running 't Hooft coupling is 



A/ (a) = gj(a)N c 



2N 2 



l3(ReTr U p /N c 



(6) 



The extra factor involving the plaquette is a standard mean-field (or tadpole) 
improved version of /3 and the naive A (a) we would derive from it. Such 
improvements are customary because the naive lattice coupling is known to 
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Figure 4. 't Hooft coupling at the scale a in physical units. 



be very poor in the sense of having very large higher order corrections. 

We extract A/ (a) and a^/a for various values of a. The latter expresses 
a in physical units so that a plot of A/ (a) against a^/a is a plot of how the 
coupling runs. If the large N c limit requires a fixed 't Hooft coupling then we 
would expect that such plots tend to a fixed curve as N c — > oo. As we see in 
Fig ^ not only does this seem to be the case, but the limit is already achieved 
at the smallest non-trivial values of N c . 



5 k-Strings 

We can consider confining strings between static colour charges in arbitrary 
representations. However gluon screening means that the effective represen- 
tation can be changed dynamically. It is therefore useful to label charges by 
their transformation properties under the centre of the group since gluons 
transform trivially under the centre. Suppose the source transforms by a fac- 
tor of z k under a global gauge transformation z belonging to the centre Zn c 
of the SU(iV c ) group. Call the lightest confining string joining sources in this 
class the fc-string. The usual string between quarks is the k = 1 string. What 
is the tension o~k of such a string as a function of k and N c ? There are some 
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Figure 5. k = 2 string tensions in SU(4) (•) and SU(5) (o). 



conjectures (see @ for details). For example, a form 



crk sm N. 



A'TT 



(7) 



has been conjectured in an M(-theory)QCE) approach to QCD U. Another 
relevant example is the old Casimir scaling El hypothesis 

cr, _ k(N c - k) 

~ ~^~T (8) 

as well. (Note we use a = u\ from now on.) 

In Fig ^ I show the lattice values of Ok=2 /cr as a function of a 2 a obtained in 
oura recent SU(4) and SU(5) calculations. (With continuum extrapolations.) 
One sees that 0^=2! & *C 2 i.e. non-trivial strongly bound fc-strings do indeed 
exist. Moreover the continuum value lies between the Casimir scaling and 
MQCD conjectures, which are numerically very similar. 

A more recent and more accurate calculation^ in SU(4) and SU(6) favours 
the MQCD conjecture. On the other hand our preliminary anisotropic lattice 
calculations □ seem to favour an intermediate value. 
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6 Corrections to the N c 



= oo Limit? 



Suppose that the fc-string tension is eventually found to satisfy Casimir scal- 
ing, as in eqn(||). Consider a fc-string wrapped around a spatial hypertorus of 
length I. This represents an energy eigenstate of the finite volume Hamilto- 
nian with mass rrik = l<Jk with corrections of 0(1/ 1 2 ) that we can neglect for 
large enough Assuming eqn(^) we obtain 

, k(N c - k) , r , k(k- 1) 
m k =lo k =la K N _ X ' =lo{k- V N ' +...}. (9) 

Consider the iV c -dependence of this, keeping I fixed in units of, say, a. We see 
that the leading correction is 0(1/N C ). This contradicts the usual diagram- 
matic result that the leading correction in SU(-/V C ) gauge theories should be 
0(1/N 2 ). 

Perhaps we should not regard a string of fixed k as being the 'same' state 
as N c is varied. Then the above would not worry us. But suppose that (part 
of) the glueball spectrum arises fmm closed strings of flux, as for example in 
the Isgur-Paton flux tube model Ej. As pointed out in til one can form such 
loops out of fc-strings as well as out of fundamental strings, leading to sectors 
of states that are scaled in mass by a k /a. Eqn(^) tells us that such glueball 
states will then have 0(1/N C ) rather than 0(1/N 2 ) corrections as N c — > oo. 

Of course we do not yet know whether fc-strings satisfy Casimir scaling 
or not. (Note that the same issue arises in D=2+l.) But the general point 
is that we have an explicit example of how non-perturbative effects - string 
formation - might lead to a violation of the usual colour counting rules. This 
is interesting whether or not reality chooses to make use of this possibility. 

7 Conclusions 

I have not had time to discuss topology. Here one finds i that the SU(oo) 
topological susceptibility is non-zero and not Jffiry different from the SU(3) 
one. This is important for our understanding E3 of how the rf gets its large 
mass. Moreover fluctuations that are unambiguously instanton-like disappear 
from the vacuum as N c grows. I have also not discussed the decaafining 
transition: the nature of this transition is being actively investigated E3i3. 

What I have shown in this talk, using fully non-perturbative calculations, 
is that the large-iV c limit is smooth, confining and is achieved precociously 
for many physical quantities. Not only is N c = 3 close to N c — oo but so is 
N c = 2. There are new stable strings at larger N c and their string tensions 
are intrigueingly close to both the MQCD and Casimir scaling conjectures. 
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One obtains the N c — > oo limit by keeping fixed the 't Hooft coupling. 
This is as expected. Not expected was the observation that the fc-strings 
provide, in principle, an explicit avenue by which states can acquire anomalous 
0(1/N C ) corrections. 

Lattice calculations which will make our knowledge of SU(A C ) gauge the- 
ories much more extensive are under way. 
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